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Abstract 

We discuss the electroweak phase-transition in the early universe, using non- 
perturbative flow equations for a computation of the free energy. For a scalar 
mass above ~ 70 GeV, high-temperature perturbation theory cannot describe 
this transition reliably. This is due to the dominance of three-dimensional 
physics at high temperatures which implies that the effective gauge coupling 
grows strong in the symmetric phase. We give an order of magnitude-estimate 
of non-perturbative effects in reasonable agreement with recent results from 
electroweak lattice simulations. 
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1. High Temperature Phase- Transitions 

There has been a lively interest in phase-transitions in gauge-theories over the last 
decades since the original work by Kirzhnits and Linde KJ, indicating that sponta- 
neously broken symmetries are restored at high temperatures. The most promi- 
nent examples are the electroweak phase-transition, i.e. the restauration of the 
SU(2)i x £/(l)y-symmetry of the standard-model, and the transition in QCD where 
(approximate) chiral symmetry is restored. These transitions are of interest for dif- 
ferent reasons. In the case of the electroweak phase-transition the most prominent 
question is whether the observed baryon-asymmetry could have been produced dur- 
ing such a transition in the early universe 0. For this scenario to work, one needs 
besides CP-violation also a sufficiently large deviation of baryon-number-violating 
processes from thermal equilibrium ||. This translates into the requirement that 
the transition be strong enough first order. On the other hand, if in the standard 
model the transition is a second or only weakly first order one, or even an analytical 
crossover instead of a true phase transition, an extension of the standard model is 
needed. Then either a B — L asymmetry could be generated in the early universe 
or the electroweak phase transition could become sufficiently strongly first order [J2J. 
An answer to the question if the observed baryon-asymmetry originated in the elec- 
troweak phase-transition requires a detailed understanding of the dynamics of the 
transition. 

The interest in the chiral transition in QCD is more directly related to speculations 
about the possibility of experimental access to the high-temperature state of QCD at 
heavy-ion colliders. At present, there is no convincing evidence that such a state has 
already been detected. A better understanding of the high temperature phase would 
certainly be helpful in answering the question of clean experimental signatures ||]. 

For QCD there was never any question about the need for non-perturbative meth- 
ods in examining the details of the transition. Correspondingly there are a host of 
studies with the help of lattice-simulations or effective models. In the case of the 
electroweak transition the hope that a description by means of high-temperature per- 
turbation theory might be sufficient prevailed for some time. For the most prominent 
qualitative features this view seems justified for a small mass of the Higgs scalar. For a 
realistic scalar mass above the experimental bound it has been argued, however, that 
a quantitative description of the high-temperature behavior and the phase-transition 
is only possible by non-perturbative methods since strong effective couplings are 
involved [5-8]. This holds despite the fact that the zero temperature electroweak 
interactions are weak. 

The deeper reason for the breakdown of perturbation theory lies in the effective 
three-dimensional character of the high-temperature field theory ||. Field theory 
at nonvanishing temperature T can be formulated in terms of an Euclidean func- 
tional integral where the "time dimension" is compactified on a torus with radius 
T _1 . For phenomena at distances larger than T _1 the Euclidean time dimension 



cannot be resolved. Integrating over modes with momenta p 2 > (27rT) 2 or, alterna- 
tively, over the higher Fourier modes on the torus (the n^O Matsubara frequencies) 
leads to "dimensional reduction" to an effective three-dimensional theory. This is 
very similar to dimensional reduction in Kaluza-Klein theories [|Hj for gravity. The 



change of the effective dimensionality for distances larger than T _1 is manifest in the 
renormalization group approach || for a computation of the temperature-dependent 
effective potential or free energy. Here one integrates over all fluctuations with mo- 
menta p 2 > k 2 and follows the dependence of the effective potential on the infrared 
scale k, finally letting k — *■ 0. The scale dependence of the effective renormalized 
couplings is governed by the usual perturbative /3-functions only for k 2 > {2tiT) 2 . In 
contrast, for k 2 < (2itT) 2 the running of the couplings was found to be determined 
by three-dimensional /3-functions instead of the perturbative four-dimensional onesQ. 
As an alternative to integrating out all modes with p 2 > (2ttT) 2 an effective three- 



dimensional theory for the long distance electroweak physics is also obtained [[L2| [13 
by integrating out the higher Matsubara frequencies^]. 

If the three-dimensional running of the couplings becomes important, the physics 
of the phase-transition is dominated by classical statistics even in case of a quantum 
field theory. A second order phase-transition is characterized by an infinite corre- 
lation length. The critical exponents which describe the behavior near the critical 
temperature are always those of the corresponding classical statistical system. Since 
the fixpoints of the three-dimensional /^-functions are very different from the four- 
dimensional (perturbative) fixpoints, we conclude that high-temperature perturba- 
tion theory is completely misleading in the vicinity of a second order phase-transition. 
This argument extends to sufficiently weak first order transitions. A second related 
example for the breakdown of perturbation theory is the symmetric phase of the 
electroweak gauge theory. The gauge bosons are massless in perturbation theory and 
the three-dimensional running always dominates at large distances |J. 

In order to understand the high-temperature behavior of a theory we should un- 
derstand the qualitative features of the /3-functions in three dimensions. These f3- 
functions have nothing to do with the ultraviolet regularization of the field theory - in 
this respect there is no difference between vanishing and nonvanishing temperature. 
They are rather related to the infrared behavior of the theory or the dependence 
of Green functions on some sort of infrared cutoff. According to Wilson's concept 
of the renormalization group these /3-functions describe the scale dependence of the 
couplings if one looks at the system on larger and larger distances. For an under- 
standing of systems with approximate scaling in a certain range it is useful to define 
dimensionless couplings. One divides out an appropriate power of the infrared cutoff 
k which plays the role of the renormalization scale. For example, the gauge coupling 
g in the effective three-dimensional theory is related to the four-dimensional coupling 



1 Related arguments in a different context can be found in pij . 

2 For an earlier treatment of dimensional reduction in high-temperature QCD see 114[ 



<74 and the temperature by 



For the S77(2)-Higgs model in three dimensions relevant for the electroweak phase 
transition, the dependence of g 2 on the scale k is givenQ by (3g2 = ^ = — -^^g\ — . . . 
with t = In k 0. One concludes that a non-abelian gauge theory like the electroweak 
theory is confining also in three dimensions. We have depicted the running of the 
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Figure 1: The running of the four-dimensional couplings g\ (solid line) and A4 
(dashed lines). The initial values of A4 correspond to scalar masses of 35 and 
80 GeV respectively. 



four-dimensional gauge coupling g\ and the four-dimensional quartic scalar coupling 
A4 in fig. 1. The deviation of gf(k) from the zero temperature value g\ = 4/9 can 
be interpreted as a measure for the validity of the one-loop approximation. For 
sufficiently small initial A4 (small physical Higgs boson mass) X^k) reaches zero for 

3 The coefficients depend on the precise choice of the infrared cutoff. 



k much larger than the three-dimensional confinement scale. One then expects a 
first-order transition which is analogous to the four- dimensional Coleman- Weinberg 
scenario [TJJ. Typical mass scales are of the order k cw where X(k cw ) = 0. In this case it 
is expected that high-temperature perturbation theory gives reliable results. On the 
other hand, if the three-dimensional confinement scale A^|j (the value of k for which 
the gauge coupling diverges or becomes very large) is reached with A(A^ nf ) > the 
behavior near the phase-transition is described by a strongly interacting electroweak 
theory. Then strong effective coupling constants appear not only in the symmetric 
phase, but also in the phase with spontaneous symmetry breaking. 

In any case, for a calculation of interesting quantities of a first order transition 
such as the critical temperature, the nucleation rate of critical bubbles, the surface 
tension, the latent heat etc. information about the strongly interacting symmetric 
phase is needed. In view of this fact, there have been different approaches to the 
electroweak phase-transition by means of lattice simulations [L7|], gap-equations 
18| , effective models and the application of renormalization group concepts || |§. 



In the following we want to review this last approach and present some comparison 
with results from 3-dimensional lattice simulations. 

2. The Average Action 

A useful tool for describing the running of couplings in arbitrary dimension is the 
average action [|I9|| . Consider a simple model with a real scalar field x- The average 
scalar field is easily defined by 

<f> k [x) = fd d yf k (x-y) X (y) (2) 



with f k decreasing rapidly for (x — y) 2 > k~ 2 and properly normalized. The average is 
taken over a volume of size ~ k~ d . The average action T k [cp] obtains then by functional 
integration of the "microscopic variables" \ with a constraint forcing <p k {x) t° equal 
the "macroscopic field" ip(x) up to small fluctuations. It is the effective action for 
averages of fields and therefore the analogue in continuous space of the block spin 
action [20] on the lattice. All modes with momenta q 2 > k 2 are effectively integrated 
out. Lowering k permits to explore the theory at longer and longer distances. The 
average action has the same symmetries as the original action. As usual it may be 
expanded in derivatives, with average potential Uk(p),p = |y? 2 , kinetic term, etc. 

r k = J d d x \u k { P ) + x -z k { P )d^d^ + ...} . (3) 

In a suitable formulation the effective average action corresponds to the gen- 
erating functional for 1PI Green functions with an infrared cutoff set by the scale k. 
It interpolates between the classical action for k — > oo and the effective action for 



k —>■ 0. In this version an exact nonperturbative evolution equation |2T] describes the 
dependence of on the infrared cutoff k (t — In k) 



dt K 2 J (2tt) 

Here Rk(q) is a suitable infrared cutoff which may depend on q 2 , as, for example, 

k 



Rk = g 2 exp (—q 2 /k 2 ) (1 — exp (—q 2 /k 2 )) 1 or R k = k 2 . The two-point function 1^ 
obtains by second functional variation of 



r ?V.9) = TT%M- (5) 

Therefore (rjj, + Rk)" 1 is the exact propagator in presence of the infrared cutoff Rk 
and the flow equation (|J) takes the form of the scale variation of a renormalization 
group- improved one- loop expression [p^j . We emphasize that the evolution equation 
is fully nonperturbative and no approximations are made. A simple proof can be 
found in |21j . The exact flow equation can be shown to be equivalent with earlier 



versions of "exact renormalization group equations" |23| and it encodes the same 



information as the Schwinger- Dyson equations 

An exact nonperturbative evolution equation is not yet sufficient for an investi- 
gation of nonperturbative problems like high-temperature field theories. It is far too 
complicated to be solved exactly. For practical use it is crucial to have a formula- 
tion that allows to find reliable nonperturbative approximative solutions. Otherwise 
speaking, one needs a description of in terms of only a few fc-dependent couplings. 
The flow equations for these couplings can then be solved numerically or by analyt- 
ical techniques. It is on the level of such truncations of the effective average action 
that suitable approximations have to be found. In this respect the formulation of the 
effective average action offers important advantages: The average action has a simple 
physical interpretation and eq. (§) is close to perturbation theory if the couplings are 
small. The formulation is in continuous space and all symmetries - including chiral 
symmetries or gauge symmetries || - can be respected. Since I\ has a representa- 
tion as a functional integral, alternative methods (different from solutions of the flow 
equations) can be used for an estimate of its form. Furthermore, the flow equation (£|) 
is directly sensitive to the relevant infrared physics since the contribution of particles 
with mass larger than k is suppressed by the propagator on the r.h.s. of eq. (Q). The 
closed form of this equation does not restrict one a priori to given expansions like 
in 1PI n-point functions. In addition the momentum integrals in eq. (H) are both 
infrared and ultraviolet convergent if a suitable cutoff Rk is chosen. Only modes in 
the vicinity of q 2 = k 2 contribute substantially. This feature is crucial for gauge the- 
ories where the formulation of a gauge-invariant ultraviolet cutoff is difficult without 
dimensional regularization. 



3. The Running Gauge Coupling 

We are now ready to discuss the running of the three-dimensional gauge coupling. 
We start from the effective average action for a pure SU(N C ) Yang-Mills theory. It 
is a gauge-invariant functional of the gauge field A and obeys the exact evolution 
equation || (with Tr including a momentum integration) 

9 t r/ll _ 1 rrx. \ dR k\ A \ M2)[.i , r gaugc(2) [ ' ^ 



dt T k [A] = -Tr | \ & [A] + r™[^] + R*[A]) | - e k [A] . (6) 

Here r| auge ^ [A] is the contribution from a generalized gauge-fixing term in a covariant 
background gauge and €k[A] is the ghost contribution M. The infrared cutoff R k is in 
general formulated in terms of covariant derivatives. We make the simple truncation 

Y k [A] = l -J ' d d xZ F>k F„ u F^ 

rr* e [A,A] = ±-jd d xZ F ^D,[A]{A»-A»)Y (7) 

with background field A and Tf nge(2) [A] = rf uge(2) [A, A = A}. In d dimensions 
the gauge coupling g appearing in F^ u and is a constant independent of k. The 
effective /c-dependent coupling can be associated with the dimensionless renormalized 
gauge coupling 

g\k) = k d -%{k) = k d ~*Z^ k g 2 . (8) 
The running of g 2 is related to the anomalous dimension t] F 

at 



dg 2 

dt ~ ^ 



(3 g 2 = (d-A)g 2 + r ]F g 2 . (9) 



Evaluating (§) for configurations with constant magnetic field and a = 1 it was found 
to obey approximately^ 



?f = (d-4)g 2 - f ^ a f 2 (10) 
dt 1 ,y l-fN c v d b d g* 

with 



v- 1 = 2 d+1 7rir 



(26-rf)(d-2) d _ 4 
a d - 44 n x 

(24-d)(d-2) d _ a 



4 In lowest order in the e-expansion |2q| the denominator in the last term is absent and 11303 is 
replaced by 114(24. 



Only the momentum integrals (x = q 2 ) 



1 2 Jo dt\dx ) 

t - -l^r^^m 1 " 1 < 12 > 

depend on the precise form of the infrared cutoff Rk appearing in 

P(x) =x + Z^Rkix) . (13) 

In four dimensions one has = 1, v 4 = 1/327T 2 and eq. (|H]) reproduces the one- loop 
result for (5 g 2 in lowest order g A . For the choice of Rk of || where 64 = 1 an expansion 
of eq. (pT0| ) in powers of g 2 also gives 93 % of the perturbative two-loop coefficient. 
We observe that the approximations leading to ([3~Q|) are valid only for \tjf\ < 1. For 
larger values of \tif\ we use a rougher estimate where bd is set to zero. 

Concerning the high-temperature field theory we should use the three-dimensional 
/3-function for k < kx, where kx = 2nT is the scale where the three-dimensional 
running sets in. The "initial value" of the gauge coupling reads (? 2 (fcr) = 20^ (/ct) 
with a w ~ ^ the four-dimensional weak fine structure constant. For k < k T the 
three-dimensional gauge coupling increases with a power behavior instead of the four- 
dimensional logarithmic behavior. The three-dimensional confinement scale A^ f - 
where g 2 diverges - is proportional to the temperature. Similarly, we may define 
(somewhat arbitrarily) the scale k np where nonperturbative effects become important 
by \r]F(k np )\ = 1. For the electroweak theory and the choices ([□]) P(x) = x + k 2 
(P(x) — x/(l — exp -p)) one finds @ 

Afl = 0.14T (0.12T) 
k np = 0.35T (0.31T) . (14) 

For the symmetric phase of the electroweak theory one therefore has to deal with 
a strongly interacting gauge theory with typical nonperturbative mass scales only 
somewhat below the temperature scale! Similar to QCD one expects that conden- 
sates like <FijF l i> play an important role [6-8]. More generally, the physics of 
the symmetric phase corresponds to a strongly coupled SU(2) Yang-Mills theory in 
three dimensions: The relevant excitations are 'W-balls" (similar to glue balls or 
strongly interacting jy-bosons) and scalar bound states. All "particles" are massive 
(except the "photon") and the relevant mass scale is set by A^ nf ~ T. For the W- 
boson masses we expect typical masses between fc np and A^. on f (cf. (|T^D). Also the 
values of all condensates are given by appropriate powers of the temperature. Since 
the temperature is the only scale available the energy density must have the same 
T-dependence as for an ideal gas 



p = cT 4 



(15) 



Only the coefficient c should be different from the value obtained by counting the 
perturbative degrees of freedomf]. We expect that quarks and leptons form SU(2) 
singlet bound states similar to the mesons in QCD^. A chiral condensate seems, 
however, unlikely in the high-temperature regime and we do not think that fermions 
play any important role for the dynamics of the electroweak phase-transition. The 
"photon" (or rather the gauge boson associated to weak hypercharge) decouples from 
the jy-balls. Its effective high temperature coupling to fermion and scalar bound 
states is renormalized to a very small value. As for the phase with spontaneous 
symmetry breaking, the fermions and "photon" can be neglected for the symmetric 
phase. We conclude that the high-temperature phase-transition of the electroweak 
theory can be described by an effective three-dimensional Yang- Mills- Higgs system. 
It is strongly interacting in the symmetric phase. Depending on the value of the mass 
of the Higgs boson it may also be strongly interacting in the phase with spontaneous 
symmetry breaking if the temperature is near the critical temperature. A more 
detailed investigation of this issue will be given in the next section. 

4. Renormalization Group improved Effective Potential for 
the Electroweak Phase- Trans it ion 

For an approximate study of the electroweak phase-transition we will now investigate 
a three-dimensional Si7(2)-Higgs model ||. It is related to the full electroweak theory 
at high temperatures by setting the [/(l)y-coupling to zero and integrating out the 
non-static modes of all fields as well as all modes of the 0-component of the gauge- 
field. For an extensive discussion of this procedure, see [13|. We will work with the 
truncation 

<P, A„ Aj = J d d x (u k {p) + Z^ k \D^\ 2 + ^Z Fjk F^F^ 

+±.Z Ftk (D,[A](A^-An) 2 ) (16) 

and a simple truncation in the ghost sector ||. Here p = (p'cp and group indices 
are omitted. The form of the average potential U k (p) is left arbitrary and has to be 
determined by solving the flow equation. For the Abelian Higgs model the evolution 
equation for the average potential was computed in the approximation ( |16|) in ref- 
erence ||. Inserting the appropriate SU(2) group factors we obtain in the Landau 
gauge (a = 0), with t = In k 

+ In (g 2 + k 2 + mfj + 3 In (q 2 + k 2 + mfj) (17) 



5 A similar remark also applies to high temperature QCD. We expect quantitative modifications 
of early cosmology due to the difference between c and the ideal gas value. 

6 We use here a language appropriate for the excitations of the three-dimensional Euclidean 
theory. Interpretation in terms of relativistic particles has to be used with care! 



where the mass terms read 



m 2 B = \z^ k glp ■ m\ = (U' k (p) + 2pU' k \p)) /Z^ k ; m\ = U' k (p)/Z^ k . (18) 

We have used here a masslike infrared cutoff R k = Z k k 2 . The partial derivative 
on the right hand side of ( |17D is meant to act only on R k and we omit contributions 
arising from the wave function renormalization Z k in R k . Primes denote derivatives 
with respect to p. 

This flow equationF] constitutes a nonlinear partial differential equation for the 
dependence of U on the two variables k and p. In our case it holds for the three 
dimensional potential U3 and a correspondingly normalized scalar field p%. They are 
related to the usual four dimensional quantities by U3 = U4/T , p 3 = p±/T. Equation 
fllTD is the basic equation of this section and has to be supplemented by corresponding 
equations for the /c-dependence of g\ and Z 9 . For given initial conditions at fcy we aim 
for a solution for k — > in order to compute the free energy U . The gauge coupling 
in ( PT5[ ) stands for the three dimensional running renormalized gauge coupling g%(k). 
Its value at the scale fc T is given by 

^Wl-f^) (19) 



where 



m 2 D = lgl(k T )T 2 (20) 



accounts for the effects of integrating out the Aq mode in lowest order ||13|| . The 
evolution equation for the running gauge coupling in the pure Yang-Mills theory has 
been given above (eq. fllOD ), and reads in lowest order 

j t ai = v# = -^tm- 1 ■ (21) 

The deviation of r from one accounts for the small contributions of scalar fluctuations 
which are not included hereQ. Equation ( pij ) is easily solved, 

1 1 23r / 1 1\ 

+ 7T7zhr-r (22) 



g!(k) g 2 3 (k T ) 24n\k T k 



7 The ultraviolet divergence on the right hand side of equation (|rj) is particular to the use of a 
masslike infrared cutoff. For d = 3 it concerns only an irrelevant constant in U and is absent for 
dU'/dt. 

8 For a suitable choice of wave function renormalization constants in the infrared cutoff for 
the gauge bosons the lowest order result becomes independent of the gauge parameter a and can 
therefore be used for the Landau gauge employed here. 



and yields the confinement scale (cf. eq.([TJ|)) in lowest order 



Ms+ssfe)" 1 - (23) 

Furthermore, we will need the anomalous dimension of the scalar field. For our 
purpose it can be approximated by |6j 

* = -nsr = -stfw*"' (24) 

and we set = 1 for k = fcy. 

A convenient quantity for an investigation of the effective potential is the in- 
dependent quartic coupling 

MP) = t£(p) = (25) 

Knowing for = the function A 3 (p) = A 3j0 (p), the high temperature effective 
potential U(p) = Uo(p) can be reconstructed by integration and translation to a four 
dimensional normalization. One of the two integration constants is irrelevant and the 
other (the mass term linear in p) can be found by adapting U(p) to the perturbative 
result for large p where the three-dimensional running of the couplings is irrelevant. 

The evolution equation for the /c-dependence of A 3 (p 3 ) can be inferred from flTTD 
by differentiating twice with respect to p 3 and reads 

dhM _ 3 / Zl k g{ k e 6Z- 2 fc A§, fc (p 3 )fc 2 2Z- 2 fc Aiy P3 )fc 2 \ 

dt 32tt^2 + m|) 3/2+ ^ + m 2 f 2 (k 2 + m 2 ) 3/2 J 1 ' 

where we have neglected terms oc U^\p 3 ) and U^(p 3 ). We report here on an ap- 
proximate solution of the flow equation for A 3 (p 3 ) for k = ||. It is based on the 
observation that the p 3 -dependent mass terms m? B) m 2 , and m 2 act in equation ( |P7D as 
independent infrared cutoffs in just the same way as k 2 . A variation of m 2 for k 2 = 
is roughly equivalent to a variation of k 2 at m 2 = 0. We use this observation to trans- 
late the flow equation (|26| ) into a renormalization group equation for A 3 (p 3 ) at k = 0: 
In equation (p6|) we replace j- by -rL = mB-^— and the factors k 2 (k 2 + m 2 ) 3 ^ 2 by 



dt "J dt' " uis dm B 

m _1 . We can then work with a new effective infrared cutoff k' = itlb which is a 
function of p 3 (we omit the prime on k in the following), 

k 2 = m 2 B = ±Z v (k)g 2 (k)p 3 . (27) 

This procedure transforms equation ( [ZED a simple differential equation for A 3 (p 3 ) = 
A 3 (&(p 3 )). In terms of the renormalized coupling 

X R (k) = Z- 2 (k)X 3 (k) (28) 



it rea 

Wt XR{k) = tkk $ {k) + + ^) ~ X R 2 ( k )9z(k)) - ^9 2 s(k)X R (k). (29) 

For the terms oc and oc ^ from equation (|26| ) we have approximated in ( jT8|) 
[/'(p 3 ) ~ p 3 A#(p 3 ) which amounts to neglecting the mass term. For negative A#(/c) 
our approximation does not describe properly the effect of the scalar fluctuations. 
Since their contribution is small in this region we simply omit the terms oc X 3 j( 2 once 
Xn(k) becomes negative. 

In order to solve the flow equation (EUp we need an initial value A 3 (£/r)> in addition 
to (0). This will depend on the (zero temperature) scalar mass M#. Furthermore, 
the p-integration of (^) leads to a term — p 2 (T)p 3 with the temperature dependent 
mass term p 2 (T) as an integration constant. The values of A 3 (£/r) and p 2 {T) describe 
how a given four-dimensional model at T > translates into an equivalent three- 
dimensional one. They can be fixed by the observation that for large p, such that 
m 2 B (p) > k^, the one-loop expression for the effective potential |26| 



UM) = -M 2 (7>3+2 (A 3 + AA 3 )p 3 2 - — (6m| + 3m|+m? + 3mi) (30) 

should be a good approximation. Here the masses are given by fh\ = 3A 3 p 3 — p 2 (T), 
and mf, = A 3 p 3 — p 2 (T), and we set Z v = 1. Also the correction 

AA - ^ (l I ^W2 Mj\ 

is chosen such that A 3 = (8n 2 T 2 / g 2 (k T )) . High temperature perturbation theory 
to one loop with the A integrated out leads then to the initial value 

We observe that the inclusion of two-loop effects or fermions will change the relation 
between A 3 (/ct) and M#. This leads to a rescaling between the scalar mass quoted 
in this work and the true physical mass. The mass term p 2 (T) has in two- loop 
perturbation theory the genuine temperature dependence 

p 2 (T) = (3M 2 H - 7 T 2 + 7T 2 In . (33) 

Here (3 and 7 are independent of T They depend, however, on the regularization 
scheme and this is particularly important for 7 as indicated already by a possible 

9 A more formal justification for eq. ( |29| ) can be found in ||. 

2 ft, \ 

10 To one loop order one has (A D integrated out, no quarks, a w = 9 ) 

1 7r ^ 3m D , Ml 

p = — , 7 = —a ' - 

and no T 2 In T term. 



change of scale in the InT-term. Some care is needed for the proper comparison 
between three- and four-dimensional lattice regularizations, high temperature per- 
turbation theory in various versions and our renormalization group approach. The 
most reliable way of comparison seems to us to equate renormalized quantities in the 
various approaches. A good candidate for fixing (3 and 7 seems to be the expectation 
value $o(T) at two different temperatures T x and T 2 

Combining equations (|25|), (27|), (|28D and fl29|) with flow equations for g\ and 
77^, we can compute the p 3 -dependence of the high temperature effective potential 
as a solution of the flow equation. It is interesting to note that except for the last 
term arising from the anomalous dimension, eq. ([29|) can also be directly obtained 
by taking appropriate derivatives of the one loop formula (|30|) , treating mass ratios 
such as mfl/mi as ^-independent and replacing at the end the couplings g\ and A3 
by running couplings evaluated at the scale k. For Z v = 1 and gf , A3 independent of 
k equation (|29|) reproduces exactly the one loop result. Our renormalization group 
improvement enables us to include the effects of running couplings and the anomalous 
dimension. This should reproduce a large part of the two-loop corrections and also 
higher contributions. The main changes as compared to the 1-loop calculations can be 
understood from the corresponding differential equations: The inclusion of 77^ lowers 
the scale at which U" changes sign, whereas the running of the gauge coupling acts 
the opposite way. Thus at a given k, corresponding to a given p 3 , the running of 
makes the potential bend up less than the 1-loop calculation predicts. The inclusion 
of the term oc A^ 2 #3 will have the same effect but is quantitatively important only 
for a large scalar mass. 

Translated to the effective potential the running of g\ should strengthen the phase- 
transition and lower the critical temperature as compared to the one-loop result. This 
is what one would naively expect, since the first order character of the transition is due 
to the gauge-boson loops. Thus, enhancing the coupling should give a transition more 
strongly first order. On the other hand, non-perturbative condensation phenomena 
may have the opposite effect. Also the inclusion of the anomalous dimension rj^ 
weakens the transition. To estimate the precise effects of the running of the couplings 
on the shape of the potential we proceed to a numerical investigation of equations 
(0), (H), and (H). (For all numerical work we use g 4 = 2/3, M w = 80.6 GeV, 
and r = 1.) Additional effects from condensation phenomena will be qualitatively 
discussed and added later. In fig.s2, 3 and 4 we show the effective potential as 
obtained with our method for different temperatures and masses of the Higgs scalar. 
In all plots we use four dimensional quantities and plot 5U = U^Q) — U^O) versus 

11 Here Ti,2 should be sufficiently below the critical temperature such that the potential minimum 
occurs in a region of p where two-loop perturbation theory is reliable. For purposes of comparison 
with three-dimensional lattice results we have chosen (3 and 7 such as to obtain the same 3>oC?i,2) 
as in two-loop perturbation theory, with T\ and T2 in a region where lattice data and perturbation 
theory are in good agreement for the prediction of the location of the potential minimum. 
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Figure 2: The effective potential for Mjj = 80 GeV with possible contributions 
from condensation phenomena. All curves are at their respective critical tem- 
perature. The curve denoted by (1) is the potential without nonperturbative 
effects (a = 0, T c = 172.8 GeV). The other curves correspond to the critical 
temperature obtained by lattice simulations Jig ], T c = 167.7 GeV, and the fol- 
lowing values of the parameters for the nonperturbative effects: 

(2) - a = 0.2928, b = 3, c = 1.1, a = 1 (solid line); 

(3) - a = 0.2885, 6 = 3, c=l, a = 2 (long dashed line); 

(4) - a = 0.305, 6 = 1, c = 1, a = 1 (short dashed line); 

(5) - a = 0.3005, 6= 1, c=l, a = 2 (dotted line). 

The lower dashed-dotted line gives the perturbative potential at the lattice- 
critical temperature. 



$ = J~pi- The first dashed-dotted curve in figure 2 (denoted by (1)) corresponds to 
the critical temperature which would be obtained from our renormalization group- 
improved approach, neglecting condensation effects. The second dashed-dotted curve 
in figure 2 gives the analogous result for a temperature corresponding to the critical 



temperature inferred from lattice simulations | I6| . We will now discuss possible 
alterations due to nonperturbative condensates and demonstrate how they could lead 
to agreement between our method and lattice results. 

The effective potential shown in figures 2 through 4 is expected to be rather 
reliable for large enough values of $ where the effective gauge coupling is not yet too 



strong. On the other hand, the truncation (TO) becomes insufficient for large g and 
we expect important modifications for $ < $ np , where (eq.s ([TJ|) and (27)) 



2k np T 



npj / 

\VF(k np )\ = 1 • (34) 

In fact, the flow equations (|10|)) (^) do not account for effects like the jy-condensation 
mentioned in the preceding section. An easy way to visualize the relevance of such ef- 
fects is the introduction of a composite 577(2)-singlet field x f° r the description of con- 
densation phenomena in the gauge sector. For example, one may choose \ cx F^F 1 ^ 
or some other (properly regularized) operator. Condensation phenomena are then 
described by the vacuum expectation value Xo- With an appropriate normalization 
of x they give a contribution to the free energy 

AU 3 = x l, AU = xlT ■ (35) 

For $ = the only available scale is k np oc T, and therefore Xo — a ^np- The 
dimensionless coefficient a is expected of order one since the gauge coupling grows 

(3) 

very rapidly for k < k np and there is not much difference between k np and A conf 
(c.f. eq. (|14D ). It is also clear that the value of xo must depend on $: For large 
$ condensation phenomena should be essentially absent since the gauge coupling 
remains small. We will parameterize the ^-dependent condensation effects by an 
additional contribution to the free energy 



AU = - 



npj \ 2 
fl'np 



3 

T (36) 



where / describes how xo depends on the effective gauge-coupling, with f(z—> oo) = 
1 and f(z) vanishing rapidly for z 1. As an example we take 

f(z) = -arctan (-b{z - c) a ) (37) 

7T \2 / 

for z > c and / = otherwise. Here c indicates for which g 2 the condensation sets in, 
and 6 is a measure how fast the condensation phenomena build up. The dependence 
of AU arises indirectly through the ^-dependence of g via the identification (|27|) . 

Condensation phenomena lower the free energy around the origin ($ = 0) and 
therefore lead to a lower critical temperature [0, |§. This is consistent with the fact 
that the critical temperature computed without AU comes out systematically higher 
than found in lattice simulations. One may even use the lattice results for T c to 
give a rough estimate of the coefficient a. Employing (ft5|),(|3"7|) with b = c = 1 we 
have adapted a such that the critical temperature coincides with the central values 



from lattice simulations. This yields a ~ 0.3 for the simulation of \IB\ at Mh = 80 
GeV. An attempt for an estimate of the size of condensation effects for Mh = 35 
GeV yields a ~ 0.38 Since the condensation effects are mainly related to the 

gauge field degrees of freedom and also the scalar contribution to the running of g 2 
is small one expects a to be independent of Mh in a first approximation ||. This 
conjecture seems to be consistent within the large uncertainties of the quoted values. 
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Figure 3: As figure 2, but with M H = 35 GeV. 

(1) - a = 0, T = 94.61 GeV (dashed-dotted line); 

(2) - a = 0.2928, b = 3, c = 1.1, a=l,T= 94.13 GeV (solid line); 

(3) - a = 0.2885, b = 3, c = 1, a = 2, T = 94.16 GeV (long dashed line); 

(4) - a = 0.305, b = 1, c = 1, a = 1, T = 94.08 GeV (short dashed line). 

(5) - a = 0.3005, b = 1, c = 1, a = 2, T = 94.10 GeV (dotted line). 

The curve denoted by (6) corresponds to the value of a as obtained from lattice 
simulations for Mh = 35 GeV jfTj: 

(6) - a = 0.381, b = 1, c = 1, a = 1, T = 93.62 GeV (solid line). 



For $ = we note that xo roughly equals the three-dimensional confinement scale 
for a ~ 0.3 — 0.4. Comparison with the curves for a = (dashed line in figure 2) 
demonstrates the importance of the nonperturbative phenomena for Mh = 80 GeV. 
Our lack of quantitative knowledge of the condensation phenomena is reflected by 
the difference of the curves for a = 1 and 2 and 6=1 and 3. For a = 2 (curves (3) 



and (5) in figure 2), AU gives no contribution to the massterm at the origin and we 
observe the minimum at $ > even for the "symmetric" phase. This would lead to 
an effective magnetic mass as described earlier in the context of gap-equations ]T8| . 

For Mh = 80 GeV we observe that condensation phenomena may weaken or 
strengthen the phase transition, depending on the values of a and b. As a general 
tendency we observe that for a fast onset of the condensation (large b) the phase 
transition becomes weaker than one would expect from perturbation theory (compare 
curves (1) and (2)). In fig. 3 we give for comparison the potential at the critical 
temperature for Mh = 35 GeV, with the same parameters a, b, c, and a as in figure 
2. In addition, we also present a curve for a = 0.381 (denoted as (6) in the figure) |I7| . 
For such a low value of the scalar mass one expects perturbation theory to be rather 
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Figure 4: As figure 2, but for Mh = 140 GeV. Here only the curves corre- 
sponding to the critical temperatures with nonvanishing condensation effects 
are given. The dashed-dotted lines give the potential at the respective temper- 
atures neglecting condensation phenomena. The parameters are: 

(1) - a = 0.2928, b = 3, c = 1.1, a=l,T= 228.5 GeV (solid line); 

(2) - a = 0.2885, b = 3, c = 1, a = 2, T = 228.5 GeV (long dashed line); 

(3) - a = 0.305, b = 1, c = 1, a = 1, T = 234.3 GeV (short dashed line). 

(4) - a = 0.3005, b = 1, c = 1, a = 2, T = 231.5 GeV (dotted line). 

For the choices of parameters corresponding to (1) and (2), the transition has 
changed to an analytical crossover. 



reliable. This is confirmed by the fact that the critical temperature T c and the 
location of the minimum $o(^c) depend only very moderately on the condensation 
effects. Only the barrier height, related to the surface tension, depends substantially 
on a, being nevertheless almost independent of the precise form of / as encoded in 
b, c, and a. For a fixed temperature the condensation effects influence the shape of 
the potential only for $ ^ 40 GeV, far away from the location of the minimum. The 
strength of the transition increases with a, independent of the shape of /. 

Finally, we also try an extrapolation to larger values of the scalar mass, as shown 
in figure 4 for Mu = 140 GeV. At the critical temperature the shape of the potential 
now depends very strongly on the shape of / and an understanding of the conden- 
sation phenomena becomes crucial even for the qualitative picture. We observe that 
all curves for a ~ 0.3 have for the symmetric phase a minimum at $o > 0. As b is 
increased, this minimum moves toward the minimum corresponding to the phase with 
spontaneous symmetry breaking. For the curve with a = 1, b = 3 the two minima 
have already melted and there remains no true phase transition (solid line). For this 
form of / one would predict an analytical crossover for Mu = 140 GeV. This illus- 
trates the speculation |J that, as a function of Mh, the critical line corresponding to 
the first order transition ends at some value Mh, c - For this critical value of the scalar 
mass the phase transition would have to be of second orderQ, with vanishing scalar 
mass at the critical temperature. By accident, this situation is realized approximately 
for a = 2, b = 3 (i.e. Mh, c ~ 140 GeV in this case). It would be very interesting to 
know the true value of the critical Higgs-mass! 

Comparing figures 2, 3, and 4 we have learned that the importance of non- 
perturbative condensation effects strongly increases with Mh- We conclude that 
perturbation theory can only give a realistic description of the phase transition for a 
small scalar mass, whereas the nonperturbative effects are of crucial importance for 
the understanding of the electroweak phase transition for realistic values of the scalar 
mass. 
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